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Abstract 

We investigate the gluon distribution functions and their contributions to the Higgs boson pro- 
duction in pp collisions in the transverse momentum dependent factorization formalism. In addition 
to the usual azimuthal symmetric transverse momentum dependent gluon distribution, we find that 
the azimuthal correlated gluon distribution also contributes to the Higgs boson production. This 
explains recent findings on the additional contribution in the transverse momentum resummation 
for the Higgs boson production as compared to that for electroweak boson production processes. We 
further examine the small-x naive /ct-factorization in the dilute region and find that the azimuthal 
correlated gluon distribution contribution is consistently taken into account, and the result agrees 
with the transverse momentum dependent factorization formalism. We comment on the possible 
breakdown of the naive /ct-factorization in the dense medium region, due to the unique behaviors 
for the gluon distributions. 



I. INTRODUCTION 



Recently, several studies have found that the transverse momentum resummation for the 
gluon-fusion processes differ from those for the quark-antiquark annihilation ( electroweak 
boson/the Drell-Yan lepton pair production) processes in the Collins-Soper-Sterman (CSS) 
framework [H, [2|, in particular, in the Higgs boson production [3] and di-photon produc- 
tion [4] processes. Similar results have been found in the context of the soft-collinear-effective 
theory formalism for the Higgs boson production jsj. These results have raised some concern 
on the derivation of the CSS formalism and the associated factorization argument for the 
gluon-gluon fusion processes. In this paper, we re-examine the transverse momentum depen- 
dent (TMD) factorization for Higgs boson production in pp collisions. We find that there is 
an additional contribution in the leading power in the TMD factorization from the azimuthal 
correlated transverse momentum dependent gluon distribution ^ . With the complete TMD 
factorization results, the CSS resummation stands for the gluon-gluon fusion process. 

Meanwhile, the azimuthal correlated gluon distribution, also referred as the linearly po- 
larized gluon distribution, has been recently discussed in the context of the transverse mo- 
mentum dependent factorization formalism in, for example, the dijet-correlation in deep 
inelastic scattering process [7], di-photon production in pp collisions jsj. This gluon dis- 
tribution will lead to the azimuthal asymmetries in these reactions, and the experimental 
observation shall provide important information on it. Moreover, in Ref. jof it was found 
that the azimuthal correlated gluon distribution has an unique behavior at small- a; from 
the saturation model [lol- 1^ calculations. This property emphasizes its special role in the 



study of the small-x gluon saturation in the high energy scattering processes, in particu- 
lar, at the planed electron-ion collider [11, In this paper, we will also examine the 



TMD factorization for the Higgs boson production at small-x, by taking into account the 
azimuthal correlated gluon distribution contribution, and compare to the result obtained 



from the well-known small-x fct-factorization formalism in the dilute region [15|, [IGJ , where 
the two gluon distribution functions are the same at small-x limit. The consistency be- 
tween these two frameworks shed important insights on the factorization property for the 
hard processes at small-x. However, in the dense medium and small transverse momen- 
tum region, the azimuthal correlated gluon distribution is different from the usual gluon 
distribution function, and we can not write the Higgs boson production as the simple naive 
/cj-factorization formalism. This indicates that the naive fcj-factorization breaks down even 
for the color-neutral particle production in the dense region in hadronic scattering pro- 
cesses as also found in Ref. [l7| , similar to the situation for the heavy quark-antiquark pair 
production process fl^ . 



The rest of the paper is organized as follows. In Sec. II, we introduce the leading or- 
der transverse momentum dependent gluon distributions, including the usual one and the 
azimuthal correlated gluon distribution function. We will also discuss the CoUins-Soper evo- 
lution for these functions, which are important for the transverse momentum resummation. 
We will show that the azimuthal correlated gluon distribution also contributes to the Higgs 
boson production in pp collisions. The CSS resummation is provided in Sec. HI, where the 
gluon distributions are calculated in terms of the integrated parton distributions at large 
transverse momentum (small b±). In Sec. IV, we extend our discussions to the small- a; region, 
where the TMD factorization and fcj-factorization formalisms are compared. We conclude 
our paper in Sec. V. 
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II. TRANSVERSE MOMENTUM DEPENDENT GLUON DISTRIBUTIONS AND 
HIGGS BOSON PRODUCTION 



The transverse momentum dependent factorization is an important step to derive the 
CSS resummation for the Higgs boson production in pp coUisions Following the Drell- 
Yan example, in our previous calculations [l9[, we have studied the factorization for the 
Higgs boson production in terms of the transverse momentum dependent gluon distributions, 
where however the azimuthal correlated gluon distribution function was not considered. In 
the following, we will find that its contribution is at the same order as the usual gluon 
distribution in the TMD factorization formalism. We consider, in general case, the Higgs 
boson production in pp collisions, 

Pa + Pb^Ho + X , (1) 

where Hq represents the scalar-Higgs boson, and the hadron A is moving +z direction and B 
along the —z. Let us first introduce the spin- independent TMD gluon distributions, which 
can be defined through the following matrix 19-22|. 



P+(27r)= 

p|F+"(r,a)4a6(r,a)/:.6c(o,o^)F:+(o)|p\ , (2) 



where F^" is the gluon field strength tensor. The light-cone components are defined as 
/c^ = (fc" ± k^) 1 1/2. In the above equation, P+ = P\ represent the light-cone momentum of 
the hadron A, and x is the longitudinal momentum fraction carried by the gluon, whereas 
/c_L is the transverse momentum. For the TMD parton distributions, the gauge link 



depends on the process [23[. In this paper, we focus on the Higgs boson production, and 

the gauge link is from the past: — oo) = Pexp {^ig d\v ■ A{\v + ^) j in the 

covariant gauge, where = A^^t^ is the gluon potential in the adjoint representation, with 
tjjj = —ifabc- In a singular gauge, a transverse gauge link at the spatial infinity has to be 
included as well. Four-velocity is an off-light-cone vector = {v~,v~^,v± = 0) with 
v~ ^ v'^ to regulate the light-cone singularity for the TMD parton distributions. With this 
parametrization, the TMD parton distributions will depend on ("^ = {2v ■ P)^/f^. (In later 
part, we also use for the TMD gluon distribution from hadron A: = {2v ■ Pa)"^ /v"^.) 
An evolution equation respect to is called the Collins-Soper evolution equation, and can 
be used to resum the large logrithms [H, llj . 

In the context of the TMD factorization, we take into account the leading power contri- 
bution in terms of P±/M where Pj_ and M are the transverse momentum and mass for the 
Higgs particle. To obtain the full result in the TMD factorization, we have to take into ac- 
count the contributions from all the leading power gluon distribution functions. The leading 
power expansion of the matrix M.^'^ in the unpolarized nucleon contains two independent 



TMD gluon distributions i, HH, HI 

1 



2k^ k" 

xg{x, kj_)g^'' + xhg{x, k^) ( \ ^ - g^^ 



(3) 



where g^ is defined as g^ = —g^'^ + {P'^P^ + P^Pj^)/PA ■ Pb- In the above parameteri- 
zations, g{x,k±) is the usual azimuthal symmetric TMD gluon distribution, and hg{x,k±) 
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the azimuthal correlated TMD gluon distribution function, hg vanishes when integrating 
over transverse momentum for the matrix A^^'^, which means there is no integrated gluon 
distribution hg{x). In the literature, this function was also called "linearly polarized" gluon 
distribution. However, in order to differ from the true linearly polarized gluon distribution 



defined for the spin-1 hadrons [2J, |25|, we prefer to use the name of "azimuthal correlated" 
gluon distribution for the spin-1/2 hadrons, following the notation of Ref. [3|. Similar func- 
tional form has been discussed in the generalized parton distribution for the gluons as well 
(see, for example, Ref. [26|). 

As mentioned above, the TMD parton distribution functions depend on the energy of the 
parenting hadrons, through the variable = (2v ■ P)^/f^ ^ 2{P^Yv~ /v^ . This equation is 
better presented in the impact parameter space [20| . 

d 

Cg^dix, b±, xC, fi, p) = {Kg + Gg)g{x, b±, xC, yU, p) , (4) 

where the gluon distribution in the impact parameter space g{x, b) is the Fourier transform 
of the TMD distribution: g{x, b±) = f cPk±e'^''^'''^g{x, k±), and Kg and Gg are soft and hard 
evolution kernels, respectively. It is straightforward to extend the above equation to that 
for the azimuthal correlated gluon distribution hg [13, 2^, 



C-^h'^g-'{x,b,xC,p,p) = {Kg + Gg)h^g''{x,b,xC,p,p) , (5) 



where hg is defined as 

/.f (x, b±) = ll d'k^e^'--"- {^-^ - g^/^ hg{x, k^) . (6) 

In particular, we find that the sum of Kg + Gg are the same for the above two equations 
and, at one-loop order, read, 

A-, + G, = -^^ln^e^-i. (7) 

TT 4 

Furthermore, the Kg and Gg obey the following renormalization group equation, 

dKg dGg ^OCsCa 

The above evolution equations will be used to perform the transverse momentum resumma- 
tion for the Higgs boson production in pp collisions. 

For the TMD gluon distribution from hadron S, we can formulate similarly. We will 
also introduce an off-light-cone vector v'^ = {v~,v~^,v± = 0) with ^ to regulate the 
associated light-cone singularity, and energy dependent variable Cl = ■ PbY /v"^. The 
same Collins-Soper evolution equations can be derived as well. 

The Higgs boson production in the gluon-fusion process can be calculated from the effec- 
tive lagrangian, 

^e// = -\9,^F^,.F^^^ , (9) 

in the heavy top quark limit, where $ is the scalar field and g^ is the effective coupling. We 
will use this effective lagrangian in the following calculations. We note the finite top quark 



4 



mass effects will not change our general discussions. The leading order perturbative calcula- 
tion produces Higgs particle with zero transverse momentum. Finite transverse momentum 
is generated by higher order gluon radiation contributions. However, at low transverse 
momentum these high order corrections introduces large logarithms of ln^(M^/Pj), which 
need to be resummed to make reliable predictions. The TMD factorization is a appropriate 
way to perform this resummation. In other words, at low transverse momentum ^ 
the transverse momentum dependence can be factorized into various leading power TMD 
parton distributions, and the higher order corrections can be factorized into the relevant 
hard factors. The hard factors in the TMD factorization does not depend on the transverse 
momentum, and the resummation can be performed by solving the CoUins-Soper evolution 
equations for the associated TMD parton distributions. The lesson from the recent stud- 
ies 0] is that, in the TMD factorization, we have to include all leading power contributions. 
In particular, the azimuthal correlated gluon distribution hg was completely ignored in the 
previous study [19]. It is straightforward to obtain this contribution in the TMD formula, 
and a similar factorization can be formulated as well. After adding this contribution, the 
Higgs boson production at low transverse momentum Pj_ <^ M can be written as, 

dMM\ -'.j d%^d%^d'e^6^'\h^ + k2± + L - p±) 



d^P^dy 



(^0 



{xig{xi, fci±)x2fi'(x2, k2±)S{i±,fi,p)H{M^,fi,p) 
X SHii±,p,p)Hh{M\ft,p)} , (10) 



where we follow the notations of Ref. [19(], ctq is the leading-order scalar-particle production 
from two gluons, ctq = vr (7^/64, and y and Pi are Higgs particle's rapidity and transverse mo- 
mentum, respectively. At low-transverse momentum, the longitudinal-momentum fractions 
Xi and X2 for the two incident gluons are related to the scalar particle's rapidity y through 
Xi = Me^ / \fS and X2 = Me~^/ v^, where S is the total center-of-mass energy squared 
S = {Pa + PbY- Ci and C2 are defined above as Ci = ^{v-Pa)/v'^ and Q = 4{v-Pb)/v'^ and p 
is a scheme- dependent parameter to separate contributions to the soft and hard factors [l9|. 
The above factorization result is accurate at leading power in P]_/M'^ at low transverse 
momentum. In particular, the interference between g and hg is power suppressed in this 
limit. 

The factorization for the first term with the usual gluon distribution follows the previous 



argument [19[, and the relevant hard and soft factors have been calculated at one- loop order. 
Similar calculations can be done for the second term in Eq. (fTOj) . In particular, at one- 
loop order, we can factorize the gluon radiation contributions to the different factors in the 
factorization formula depending on the kinematic regions of the radiated gluon. For example, 
if the radiated gluon is parallel to the incoming hadron A, we factorize its contribution to 
the TMD gluon distribution hg from A. If it is parallel to the hadron B, we include that 
contribution to the TMD gluon distribution hg from B. When the gluon momentum is soft, 
we factorize its contribution to the soft factor. The hard factor is calculated from the hard 
gluon radiation in the virtual diagrams, because the real gluon radiation is power suppressed 
if all momentum components are hard in order of M. We have done the explicit calculations 
to show this factorization at one-loop order 2^ . 



It is easy to find that the soft factors for the above two terms are identical: Sh{i±) = 
S{i±). This is because the soft gluon radiations do not depend on the spin/polarization, 
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and are defined identically for these two terms. In the one-loop calculations, the hard 
factor are extracted from the virtual diagrams for the cross section and parton distribution 
calculations, and we find that the hard factors are also the same, Hh{M'^) = H{M'^) (29| . 

i^«(M^/.^p) = /7«(M^;.^p) 

O^sCa 



. f^n 1 1 2 3\ 3 , 2 1 1 2 2 2 ^ 



where a special coordinate system has been chosen in which xf^f = x\C,2 = pM"^. In a recent 
calculation for the single-spin dependent observable, a similar conclusion was also found [Hj, 
which may indicate that all the hard factors in the TMD factorization are independent of 
the spin/polarization. 

It is convenient to write down the TMD factorization formula in the impact parameter 
space. 



^P,,ly ^«Jj^^-'-'^W(...-.AM^), (12) 

where W contains contribution from the two terms in the TMD factorization, 

W{xi,X2,b,M^) = Wg{x,,X2,b,M^) + WH{xi,X2,b,M^) , (13) 

and Wg and Wh represent the contributions from the usual gluon distribution g{x,b) and 
the azimuthal correlated gluon distribution hg{x,b), respectively, 

Wg{xi, X2, b, Q^) = S{b, fi, p)H{M^, p, p)xig{xi, b, p, pM'^, p)x2g{x2, b, p, pM'^, p) (14) 
Wh{x^, X2, 6, Q') = 25(6, /i, p)H{M\ /i, p)xi/if (a;i, 6, /i, pM\ p)a;2/if (0:2, 6, p, pM\ p)(^5) 

where the 1^/1 comes from the specific tensor structure in the factorization formula Eq. (ITOl) . 
The convolutions in the transverse-momentum space now reduce to products in the impact 
parameter 6-space. In the factorization formula, the large logarithms will show up as In M^fe^ 
in the various factors in the above equations. We need to resum these large logarithms. 



III. RESUMMATION 

The large logarithms in the factorization formulas in the last section are resummed by 
following the Collins-Soper-Sterman method. The two terms in Eq. f|T3l) satisfy the CoUins- 
Soper evolution equation separately. 



dWg^r.{x^AM^) 
(91nM2 



{K + G')Wg^h{x,AM'') , (16) 



where K and G' are soft and hard evolution kernels. Since the two gluon distributions obey 
the same Collins-Soper evolution equation and the hard factors are the same, the evolution 
kernels are the same as well. Combining the Collins-Soper evolution equations for the TMD 
gluon distributions of Eqs. fl5|7|) and the hard factors at one- loop order of Eq.([TT]), we find 
that, 

ir + G' = -^ln(^e^--^o^ , (17) 
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where the p dependence between various terms cancels out. Solving the above evolution 
equations, we obtain, [l| 

6, M^) = e-^i-(^^''''^^/^^)iy,,,(x„ 6, Cl/Cl/h') , (18) 
where the large logarithms are included in the Sudakov form factors. 



5, 



Sud 



Cf/b2 /^^ 



ln(^^)A(Ci,^) + i?(Ci,C2,/i) 



(19) 



Here Ci and C2 are two parameters of order one. The functions A and B can be expanded 

00 . 00 . 

perturbatively a^, A = A^'^^ (^) and B = B^^^ (^) . Because the A coefficients come 

1=1 ^ 1=1 ^ 

from soft factor which are the same for the two terms Wg and Wh, we expect A will be 

the same as well. On the other hand, B coefficients come from the hard factors in the 

TMD factorization formulas. Therefore, they could be different 30|]. However, our one-loop 

calculations lead to the same hard factors and the same B coefficients for Wg and Wh- We 

expect that the effects discussed in Ref. js^] do not affect our calculations, and we conjecture 

that the B coefficients will be the same for these two terms at higher orders too. With this, 

we can combine the above two terms together as, 

W{x„ b, M') = e-^su,(M\i>,c,/c,) ^^^(^^^ Cl/Cl/h^) + Wh{x,, b, Cl/C'Jb')] , (20) 

where Ssud represents the universal Sudakov form factor for the Higgs boson production. 
Up to the one-loop order, we have verified this result. 

The last step of the complete CSS resummation is to formulate the Wg and Wh of the right 
hand side of Eq. (12U]) at lower scale Cf/C|6^ in terms of the integrated parton distributions. 



Wg{x.,b,Cl/Cl/b') = J2 [ ^^x[Mx[,p)x'J,{x'„p) 



^3 



xQ,(xi/x;, CJC2/b/p)Cg,j{x2/x'^, CJC^/b/p) , (21) 
Wh{x,,b,CllCllb') = Y, I ^^x;/,(a:;,/i)4/,(4,/i) 

J Xi X2 

''J 

xCh/iix,/x[, C,/C2/b/fi)Ch/jix2/x',, C,/C2/b/fi) , (22) 

where fij reprent the integrated quark/gluon distribution functions. For integrated gluon 
distribution, there is only the usual one, whereas the counterpart of hg does not exist. The 

C = J2 C^'^ (v)' coefficient funct ions can be calculated perturbatively. The coefficients Cg/i 

have been calculated up to two-loop order [3|, where Ch/i are also calculated up to order. 
In their calculations, the cross section W{b) = Wg{b) + Wh{b) in the impact parameter {b±) 
space is written as perturbative expansion of Og, from which the relevant coefficients are 
extracted by comparing with the Eqs. f l?ip2]) . In the following, we will show how we can 
calculate Ch/i from the TMD factorization formula Eq. ( !T5|) . 

To calculate the Wh in Eq. (!22|) . we compute the azimuthal correlated gluon distribution 
h'^'^ in terms of the integrated quark/gluon distribution functions and substitute into the 
factorization formula Eq. (fT5|) . First, we write down a similar factorization form for h'^'^{b±), 

1 / 9h^ \ r fix' ~ 
K^ixM) = 2 hT - j / -^Ch/.{x/x',b^,p)x'Mx\p) , (23) 
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where the pre-factor of | (^g'^ — -^r^^ comes from the basic Lorentz structure for this 

functioE0- We know that there is no integrated hg gluon distribution, which immediately 
leads to the zeroth order of expansion of the above equation vanishes. As a consequence, 
the zeroth order of Ch/i in Eq. fl22|) vanish as well, 

C-S", = = . (24) 

At order of a^, we can generate the azimuthal correlated gluon distribution from the in- 
tegrated quark/gluon distribution functions. For example, the contribution from the inte- 
grated gluon distribution is, 

hg{x,k^) = ^Ca^ J ^^-^9{x') , (25) 

where ^ = x/x'. The Fourier transform into the impact parameter space leads to. 

This Fourier transform does not generate any divergence, which is consistent with the fac- 
torization formula of Eq. f l25]) . Because the non-zero leading order expansion of Eq. fl25]) 
is at order a^, the right hand side is associated with the leading order gluon distribution, 
and there is no collinear divergence. An interesting consequence is that the non-zero leading 
order coefficients do not depend on the factorization scale [3|]. However, from the factoriza- 
tion formula Eq. ( 123|) . at order of a^, we will find out the Fourier transform will lead to a 
collinear divergence which shall be absorbed into order splitting of the integrated gluon 
distribution function. This indicates that order coefficients Cj^^. (and consequently the 

(2) 

following C^/j) will depend on the factorization scale. 

Similar results are obtained for the azimuthal correlated gluon distribution in terms of 
the integrated quark distribution with the color-factor Cp instead of Ca, 

/C(-.*)4(.r-^)>/|^i^,M. (27) 

Combining Eqs. f l26|27|) with Eq. f lT5|) . we obtain, 

CH/, = CFil-0, Cil = CA{l-0, (28) 

which reproduces the relevant resummation formula in Ref. 

From the above derivation, we find that the different resummation formalism for the 
gluon-gluon fusion processes as compared to that for the Drell-Yan lepton pair production 
process comes from the fact that there are two independent TMD gluon distribution func- 
tions at the leading order which contribute to the Higgs boson production at the same order 
in the limit of Pi ^ M. Although there are perturbative at different order in terms of the 
integrated quark/gluon distribution functions, we have to take into account the contribu- 
tions from both functions in order to completely describe the Higgs boson production at low 
transverse momentum ^ M. In particular, in certain kinematic region such as small-x 
region discussed in next section, the azimuthal correlated gluon distribution is as important 
as the usual one, where we have to include its contribution. 



^ In order to keep this factor traceless in the dimensional regulation calculations, the factor 2 in the bracket 
should be replaced by c? — 2 where d = 4 — e denotes the dimension. In our following calculations of 
Eqs. (|26l27p . since there are no divergence in the Fourier transform, we use d — 4. 
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IV. SMALL-x /ct-FACTORZATION 



The two TMD gluon distribution functions at small-x have unique properties as recently 
discussed in Ref. |9|. Applying these results into the factorization formula of Eq. ( ITU]) , we 
will be able to study the factorization of Higgs boson production in small- a; region, where 



the well-known /crfactorization formalism [10| has been applied too. In this section, we 



discuss the Higgs boson production in the small- a; region, and examine the so-called naive 
/ct-factorization approach in this kinematic region. 

From Eqs. (12^ . we notice that the gluon splitting contribution to the azimuthal corre- 
lated gluon distribution has the same small- a; enhancement as the usual gluon distribution. 
Therefore, we expect the similar BFKL evolution for hg{x,k±) at small-x in the dilute 



exp 

regime [l^, l3ll-33|. since the operator definition of hg{x, k_i_) at low- a; is also related to the 
quadrupole. As a consequence, the azimuthal correlated gluon distribution will be as im- 
portant as the azimuthal symmetric one in this kinematic limit. In particular, from the 
saturation model calculations of Ref. [9|, we know that the azimuthal correlated gluon dis- 
tribution function is the same as the usual gluon distribution function at small-x in the dilute 
region with kj_ ^ Qg, where Qg is the characteristic scale in the saturation model. This is 



also consistent with the expectation from the BFKL evolution for these two functions [14 
Therefore, in this region, we can combine the two contributions in the factorization formula 
Eq. (fTOj) into one. 



ay a ij_ J "^i±'^2± 

where we have used g{x^ k±) to represent both g and hg distribution functions, and neglected 
higher order corrections from the hard and soft factors in Eg. fflQl). The above result is 
exactly the same as that obtained in the naive-fcf factorization |16| by taking the small 



transverse momentum limit Pi <^ mH. By using the proper physical gluon polarization (k 
one automatically takes into account the contribution from the azimuthal correlated gluon 
distribution in the naive fc^-factorization approach. 

However, in the dense medium (large nucleus or extremely small-x), and in particular 
when k± ~ Q^, the azimuthal correlated gluon distribution is different from the usual gluon 
distribution [9] if they follow the definitions in Eq. (^!^. They are appropriate definitions 
for the gluon distribution in the Higgs boson production process Therefore, we can 



not combine these two terms into one universal structure as suggested in the naive /cr- 
factorization at small-x. This indicates that the naive-kt factorization breaks down even 
for the color-neutral particle production in the dense medium in the hadronic scattering 
processes. Similar conclusion has also been drawn for the rj' particle production in pA 



collisions in the saturation model calculations [17(|. However, because of the large Higgs 
mass, we should be able to modify the naive-fct factorization to establish an effective kf- 
factorization for its production at low transverse momentum Pj^ <S M, following the similar 



study in Ref. [23|. This will lead to consistent results as the TMD factorization of Eq. f llOp 
with the small-x gluon distributions calculated in the dense region. An explicit calculation, 
including high order corrections, will be very important to investigate the QCD factorization 
property for the hard processes at small-x. 



^ The difference between the resuhs in Refs. 15, 1^ vanishes in this hmit. 
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V. CONCLUSIONS 



In summary, in this paper, we have investigated the transverse momentum dependent 
gluon distribution functions and the Higgs boson production in pp coUisions in the transverse 
momentum dependent factorization approach. We found that the azimuthal correlated gluon 
distribution contributes to the Higgs boson production in the leading power of Pt/M. After 
taking into account this contribution, we will be able to explain recent findings on the 
resummation for the Higgs boson production at moderate transverse momentum. It will 
be interesting to extend this study to the di-photon production process and the associated 
resummation formalism 0, Isf . 

We further extended our discussion to the small-x region, where we compared the TMD 
factorization result with the well-known naive fcj-factorization result, and found that they 
are consistent in the dilute region. We expect they will differ in the dense region, which 
may indicate the naive /cj-factorization is violated even for the neutral particle production 
at small-x region. 
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